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In general, the stability of a band crossing point indicates the presence of a quantized topological 
number associated with it. In particular, the recent discovery of three-dimensional Dirac semimetals 
in NaaBi and CdaAs 2 demonstrates that a Dirac point with four-fold degeneracy can be stable as long 
as certain crystalline symmetries are supplemented in addition to the time-reversal and inversion 
symmetries. However, the topological charges associated with NasBi and CdaAs 2 are not clarified 
yet. In this work, we identify the topological charge of three-dimensional Dirac points. It is found 
that although the simultaneous presence of the time-reversal and inversion symmetries forces the 
net chiral charge to vanish, a Dirac point can carry another quantized topological charge when an 
additional rotation symmetry is considered. Two different classes of Dirac semimetals are identified 
depending on the nature of the rotation symmetries. First, the conventional symmorphic rotational 
symmetry which commutes with the inversion gives rise to the class I Dirac semimetals having a pair 
of Dirac points on the rotation axes. Since the topological charges of each pair of Dirac points have 
the opposite sign, a pair-creation or a pair-annihilation is required to change the number of Dirac 
points in the momentum space. On the other hand, the class II Dirac semimetals possess a single 
isolated Dirac point at a time-reversal invariant momentum, which is protected by a screw rotation. 
The non-symmorphic nature of screw rotations allows the anti-commutation relation between the 
rotation and inversion symmetries, which enables to circumvent the doubling of the number of Dirac 
points and create a single Dirac point at the Brillouin zone boundary. 


I. INTRODUCTION 

After the discovery of graphene, a class of materi¬ 
als, dubbed Dirac semimetals, have come to the fore 
of condensed matter research. In general, a Dirac 
semimetal has several Fermi points around which pseudo- 
relativistic linear dispersion relation is realized. This 
pseudo-relativistic energy dispersion forces the density 
of states on the Fermi level to vanish without opening 
of an energy gap, which is the unique property of Dirac 
semimetals distinct from ordinary metals or insulators^. 
In particular, the recent theoretical prediction^ and the 
experimental confirmation^— of three-dimensional (3D) 
Dirac semimetals in NaaBi and CdaAs 2 demonstrate that 
there are a variety of materials realizing Dirac semimetals 
in both two dimensions and three dimensions. Such a di¬ 
versity of Dirac materials requires us to find a systematic 
way to characterize and classify them. 

In general, the stability of nodal points in a Dirac 
semimetal has topological origin. This is because there 
is no characteristic energy scale, such as the Fermi en¬ 
ergy or the energy gap, characterizing the perturbative 
stability of the system. For instance, a nodal point in 
graphene carries a quantized pseudo-spin winding num¬ 
ber that is defined on a loop encircling the Dirac pointi 
On the other hand, a nodal point in 3D Weyl semimet¬ 
als is endowed with a Chern number defined on a two- 
dimensional (2D) closed surface surrounding the Weyl 
point— The presence of such a quantized topological 
charge carried by a nodal point guarantees its stability, 
hence a nodal point can be annihilated only by collid¬ 


ing with another nodal point with the opposite topo¬ 
logical charge as long as the symmetry of the system is 
preserved. Recently, there have been several theoretical 
studies which attempt to classify the topological invari¬ 
ants of nodal points, and to extend the concept of topo¬ 
logical band theory to gapless systems such as semimetals 
and nodal superconductors^—. However, in our opin¬ 
ion, a proper definition of the topological charge of Dirac 
points in NasBi or CdaAs 2 has not been given so far. 

Dirac points in NasBi or CdsAs 2 are protected by 
the time-reversal (T), the inversion (P), and the rota¬ 
tion symmetries^iii— The experimental observation of 
Dirac points in these systems demonstrates their stabil¬ 
ity, hence the presence of topological invariants associ¬ 
ated with them. Moreover, the theoretical observation of 
pair-annihilation and pair-creation of Dirac point s 22 ’ 33 
indicates that the topological charges of the two Dirac 
points should have the opposite sign. Then the ques¬ 
tion is what the nature of the topological charge is as¬ 
sociated with the Dirac points. Considering the two- 
dimensionality of the sphere surrounding a Dirac point, 
the natural candidate is either a Chern number similar 
to the case of Weyl semimetals, or a Z 2 invariant associ¬ 
ated with T symmetry satisfying T 2 = —1. However, the 
simultaneous presence of the time-reversal and the inver¬ 
sion symmetries forces the Berry curvature to be zero at 
each momentum, hence the Chern number of the Dirac 
point, which is basically the integral of the Berry curva¬ 
ture, also vanishes. Moreover, a Dirac point can carry 
a Z 2 topological charge only in the presence of SU(2) 
spin rotation symmetry together with time-reversal and 
inversion symmetries satisfying (TP) 2 = 1 as shown in 
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Ref. [ 23 ). These indicate that a special care is required 
to find the topological charge of a Dirac point in Na 3 Bi 
or Cd 3 As 2 , which should obviously be distinct from the 
monopole charge of Weyl semimetals. 

In this work, we will show that the Dirac points in 
Na 3 Bi or Cd 3 As 2 are characterized by topological in¬ 
variants of zero-dimensional subsystems defined on the 
rotation axis. Since the rotation eigenvalue is a good 
quantum number on the rotation axis, a zero-dimensional 
topological invariant can be defined by comparing the ro¬ 
tation eigenvalues of the valence and conduction bands 
at two points enclosing a Dirac point. We find that the 
nature of 3D Dirac semimetals strongly depends on the 
nature of the rotation symmetry. Namely, the ordinary 
symmorphic rotation symmetry commuting with the in¬ 
version symmetry always creates a pair of Dirac points 
having the opposite topological charges, and generates 
class I Dirac semimetals. Both Na 3 Bi and Cd 3 As 2 be¬ 
long to this class. 

On the other hand, we find a different class of Dirac 
semimetals when the system has a screw rotation sym¬ 
metry. In general, the presence of non-symmorphic sym¬ 
metries, such as screw rotations and glide mirror sym¬ 
metries, guarantees a nontrivial band connection at the 
Brillouin zone boundary^—. Also, it is proposed that 
when the double space group of non-symmorphic crystals 
satisfies certain conditions, a Dirac point can be realized 
at the Brillouin zone boundary^. Consistent with these 
results, our theoretical study shows that when the band 
degeneracy at the zone boundary is compatible with the 
time-reversal and inversion symmetries, a single isolated 
Dirac point can be created on the rotation axis. Based 
on this observation, we define a class II Dirac semimetal 
which is protected by a screw rotation symmetry and the 
inversion, which are mutually anti-commuting, in addi¬ 
tion to the time-reversal symmetry. The partial transla¬ 
tion associated with a screw rotation adds a U(l) phase 
to the rotation eigenvalue, which varies on the rotation 
axis. This projective nature of a screw rotation enables 
to circumvent the doubling of Dirac points, and create 
a single isolated Dirac point at a time-reversal invariant 
momentum at the Brillouin zone boundary. 

The rest of the paper is organized in the following 
way. We describe the general idea to define a topologi¬ 
cal charge in systems with rotation symmetry in Sec. [TT] 
Based on this general idea, class I Dirac semimetals 
are defined and systematically classified in Sec. IIIII In 
particular, we show that the doubling of Dirac points 
is unavoidable in class I Dirac semimetals. Sec. hyi is 
about the nontrivial band connection generated by non- 
symmorphic screw rotation symmetries. Here we show 
that the partial translation associated with a screw rota¬ 
tion induces a momentum dependent U(l) phase factor 
to the rotation eigenvalue, which enables to circumvent 
the fermion number doubling and protects a single Dirac 
point at the Brillouin zone boundary. Based on the dis¬ 
cussion in Sec. El class II Dirac semimetals are defined 
and systematically classified in Sec. El We present the 


conclusion and discussion in Sec. eh In the Appendix, 
we prove that there is no stable Dirac semimetal in sys¬ 
tems only with the time-reversal and inversion symme¬ 
tries based on K theory approach. The classification of 
Dirac semimetals in C 2 invariant systems shown in the 
main text is also confirmed by using the K theory ap¬ 
proach. Finally, we present a short discussion about the 
stability of 2D Dirac semimetals protected by two-fold 
screw rotations. 


II. GENERAL IDEA: ROLE OF ROTATIONAL 
SYMMETRY IN SYMMORPHIC CRYSTALS 

In general, electronic systems having only the time- 
reversal ( T ) and inversion (P) symmetries cannot sup¬ 
port a stable Dirac point with a quantized topological 
charg o 11 ! 34 : 35 . In Appendix A, we have revisited this 
known fact in a different perspective and proved it by us¬ 
ing K theory approach. Thus additional crystalline sym¬ 
metries play a crucial role to stabilize Dirac senrimetals 
realized in Na 3 Bi and Cd 3 As 2 . Here we consider the role 
of the additional rotation symmetry (Cn) in addition to 
P and T. For convenience, we first focus on 3D crystals 
with a symmorphic space group symmetry in which the 
point group can be completely separable from pure trans¬ 
lation operations. Also we choose the z axis as the axis 
for Cn rotation with N indicating the discrete rotation 
angle of 2n/N ( N = 2,3,4, 6). Under the operation of 
the Cn symmetry, the Hamiltonian satisfies 

C N H(k x ,k y ,k z ){C N )~ 1 = H(k x ,k y ,k z ), ( 1 ) 

where ( k x , k y ) is obtained from 2tt/N rotation of ( k x , k y ), 
i.e., ( k x + iky) = ( k x + ik y )e 2m ^ N . The symmetry oper¬ 
ators satisfy 

T 2 = -1, (C N ) N = -1, P 2 = 1, (2) 

and 

[T,P) = 0, [T,CV]=0, [P, Cn] = 0, (3) 

where we have considered the fact that an electron is a 
spin-1/2 particle. 

Now let us explain the general idea of how to deter¬ 
mine the topological charge of a Dirac point locating 
at a generic point fco = (0, 0, fc°) on the rotation axis. 
To determine the topological charge, we first consider a 
sphere in the momentum space surrounding the Dirac 
point at fc = fco in a Cn symmetric way as shown in 
Fig. [Q Namely, the center of the sphere sits on the rota¬ 
tion axis. At every point on the sphere, the Hamiltonian 
is invariant under the compound antiunitary symmetry 
PT satisfying (PT) 2 = —1. Moreover, the intersection 
of the k z axis and the sphere consists of a north pole 
at k N = (0,0, k^) and a south pole at k s = (0,0,/cf), 
which are invariant under the rotation. We define the 
topological charge of the Dirac point from the topologi¬ 
cal numbers associated with these two points. Since the 
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FIG. 1. Local geometry around a Dirac point at k = fco 
sitting on the rotation axis. The sphere surrounds a Dirac 
point at the center. k N and k 3 mark the points on the sphere 
crossing the rotation axis. 

Hamiltonian commutes with the rotation operator Cn at 
these points, 

[H(k N ' s ),C N ] = 0, (4) 

hence H(k N / s ) can be block-diagonalized in the 
eigenspace of Cn with the eigenvalues J m given by 

J m = exp ^ i 1 7 r^ , m = 0 ,...,JV-l. (5) 

Since the PT symmetry is not satisfied in each Cn 
eigenspace with a given J m in general, (one exceptional 
case is shown in Sec. IIII Bl) . each diagonal block of 
H(k N ' S ) belongs to the symmetry class A in terms of 
the Altland-Zirnbauer classification scheme^, hence car¬ 
ries an integer topological number Here n^ or 

indicates the topological invariant of a zero-dimensional 
system belonging to the symmetry class A, which is de¬ 
fined as 

nZ=l[N c (J m ,k N )-N v (J m ,k N )], (6a) 
n s m = ^[N c (J m ,k s )~N v (J m ,k s )], (6b) 

where N c / v (J mi k) denotes the number of conduction 
bands (c) or valence bands (v) with the eigenvalue J m 
at the momentum k. It is worth to note that a triv¬ 
ial conduction or valence band with a constant energy 
can always be added to each J m sector, so that the sum 
n m + n-m can be changed freely. Therefore the nontrivial 
topological number in the J m sector is determined by the 
difference 

"m = n* - n^. (7) 


Next let us consider possible constraints to the al¬ 
lowed v m values. For a gapped system defined on the 
sphere, the number of conduction bands and that of va¬ 
lence bands are constants independent of the momentum 
on the sphere, which leads to the following constraint 

m m 

= {[Nc(Jm, k N ) - N c (J m , k s )] 

m 

- [N v (J m ,k N ) - N v (J m ,k s )}} 

= 0 . ( 8 ) 

Moreover, since PT symmetry imposes additional con¬ 
straints between different v m values, the number of in¬ 
dependent topological invariants depends on the details 
of the symmetry as shown in Sec. nm and El However, 
as long as a Dirac point possesses a nonzero v m value, 
it guarantees the stability of the relevant Dirac point. 
Hence the set of nonzero v m can be considered as a topo¬ 
logical invariant characterizing a stable Dirac point. 

III. CLASS I DIRAC SEMIMETALS 

Class I Dirac semimetals are protected by the ordinary 
rotation symmetry commuting with an inversion symme¬ 
try, i.e., 

[P, C N ] = 0. (9) 

Let us consider the eigenstate | ip m ) of the Cn operator 
with the eigenvalue J m . The PT symmetry requires 

C N PT\^ m ) = PTC N \lpm) 

= PT J m \ij) m ) 

= r m PT\ii> m ) 

= JN-m-lPT\ll ) m ), ( 10 ) 

from which we find PT\ip m ) is an eigenstate of Cn with 
the eigenvalue Jjv-m-i- Therefore when a state with the 
eigenvalue J m is occupied (or unoccupied), there should 
be another occupied (or unoccupied) state with the eigen¬ 
value Jn-tu- i, which leads to the constraint 

Vm — VN-m — 1- ( 11 } 

For further analysis, we distinguish two cases based on 
the parity of N as shown below. 

A. Cn symmetric systems with even N 

Due to the PT symmetry, the Cn eigenspaces with 
the eigenvalues J m and Jyr-m-i can be paired as 
{Jm,JN-m- 1 } with m = 0,..., N/2- 1. [See Fig. [2] (a).] 
Since PT interchanges eigenspaces within each pair, PT 
is not a symmetry in each eigenspace separately, and each 
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(a) N = even (b) N = odd 



FIG. 2. Constraints on the rotation eigenvalues (a) for even 
N, (b) for odd N. The (black) solid circle indicates a unit 
circle in the complex plane, and each (red) dot on the circle 
denotes J m . Two dots connected by a dotted line are related 
by the PT symmetry, hence only one of them is independent. 


eigenspace belongs to the symmetry class A. Therefore an 
integer topological invariant defined in Eq. (Rral) can be 
computed in each eigenspace with J m . Considering the 
constraints shown in Eq. (jSJ) and Eq. m, we conclude 
that the topological charge of the Dirac point is given by 

(TO) • • ■ j vn__ 2 ) £ Z 2 1 . (12) 

Therefore, the topological charge of a Dirac point with 
C 4 or Cq symmetry is an element of Z or Z 2 , respectively. 
At the same time, it implies that a C 2 invariant system 
cannot support a stable Dirac point. 


B. Cn symmetric systems with odd N 

When N is odd, the PT symmetry pairs the 
eigenspaces of CV in a slightly different way as compared 
to even N cases as shown in Fig. [2] (b). At first, we find 
l(A r — 1 ) pairs of eigenspaces 

N - 3 

{Jmi Jn— m— 1} 1 m — 0, . . . , “ , (13) 

On the other hand, the remaining eigenspace with the 
eigenvalue J(aj_i)/2 is invariant under the PT symmetry, 
hence belongs to class AIL Thus, in a block-diagonalized 
Hamiltonian H(k N / s ), there are 1 (N — 1) blocks be¬ 
longing to class A and an extra block with the eigen¬ 
value J(jv_ i)/2 belonging to class AIL In both symmetry 
classes, zero dimensional systems have an integer topo¬ 
logical invariant, which is defined as the difference in the 
number of conduction bands and that of valence bands. 
Hence as in the even N case, the topological charge can 
be defined as 

Vm = n I Z l - n^, (14) 


in each eigenspace with the eigenvalue J m . The con¬ 
straints to the topological numbers v m are 

N - 3 

{i)v m = VN-m- 1 , m= 0,...,—-—, (15) 

(ii) ^2 V rn = 0. (16) 

m 

Thus the independent topological charge for a Dirac 
point is given by 

(dj, • ■ ■, v (n-z)/i) £ ^ 2 ■ (17) 

Therefore a Dirac point with C 3 symmetry has an integer 
(Z) topological charge. 


Cn 

Topological charge 

c 2 

Not allowed 

c 3 

Z 

c 4 

Z 

C 6 

Z x Z 


TABLE I. Summary of topological charges of class I Dirac 
semimetals. 


C. Applications: classification of stable Dirac 
points in 4-band systems 

Let us apply the general theory developed above to 
minimal 4-band models, and classify stable Dirac points. 
In a 4-band model, a pair of doubly degenerate bands 
cross at a Dirac point which we assume to sit on the Fermi 
level. On the rotation axis fc = (0,0, k z ), each band is 
assigned with a quantum number J m . Since the pair of 
degenerate bands should have different rotation eigenval¬ 
ues to generate a stable Dirac point, each band with the 
rotation eigenvalue J m satisfies n ^ = —n^- Namely, a 
band which is below (above) the Fermi level at the mo¬ 
mentum k N should be above (below) the Fermi level at 
the momentum k s to have a Dirac point in between. 

1. C 2 symmetric systems 

J m - 0 = exp(i 7 r/ 2 ) = i and J m = 1 = exp( 3 * 7 r/ 2 ) = 
—i = J m --\ are the only allowed C 2 eigenval¬ 
ues. Due to the PT symmetry, a pair of eigenstates 
{|V'm=o(fc))i |f/ , m=i(fe))} are always degenerate locally at 
each momentum fc, hence n ^’-fo = n m =1 an( i A m =o = 
N m —i. Then a 4-band model can be constructed by intro¬ 
ducing two pairs of eigenstates {|^_ 0 (fc)), |^_i(fc))} 
and {|'i/’m = 0 (fc)), |V’m= 1 (^))}, where A,B indicates the 
valence band (v) or the conduction band (c), respectively. 
It is straightforward to show that n^ s = 0 (m = 0, 
1 ), because if one state is occupied, among {|^(fc)), 
iV’m(Jc))}, the other state is unoccupied. Therefore 
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v m=o,i = 0 and there is no stable Dirac point with a 
nontrivial topological invariant in systems with C 2 sym¬ 
metry. 

2. C 3 symmetric systems 

Possible C 3 eigenvalues are J m =o = exp(*7r/3), 

J m = 1 = exp(z7r), J m -2 = exp(5z7r/3). Due 

to the PT symmetry, {\ip m =o{k)),\4> m= 2 (k))} and 
{\ip m= i(k)), |?/> m=1 (A;)}} form degenerate pairs, where 
I’ipm^ik)} = PT\i/j m —i(k)). Thus v m = 0 = 2 . Since 

0 + ^ 777=1 + ^m =2 = 0, we have 

^1 = ^ 777=0 = ^ 777=2 = ~^m= 1 - ( 18 ) 

Hence there is only one independent topological num¬ 
ber, Vi £ Z. Since the topological charge of a 
Dirac point can be nonzero only when the valence and 
conduction bands have different rotation eigenvalues, 
a 4-band model can be constructed by using a ba¬ 
sis Wm=o( fc ))> IV 4 = 2 (* 0 }: l^m=i( fc )>} where 

A = v (c) and B = c (v). Since = 

~\ n m =1 = ±5 and for 4-band models, the 

Dirac point has a nonzero topological invariant 

z7i = ±l. (19) 


3. C 4 symmetric systems 

Possible C 4 eigenvalues are J m =o = exp(* 7 r/ 4 ), J m -i = 
exp( 3 i 7 r/ 4 ), J m=2 = exp( 5 i 7 r/ 4 ), J m=3 = exp( 7 i 7 r/ 4 ). 
Due to the PT symmetry, {\4> m =o(k)), \tp m= 3 (k))} and 
{\ 4 ’m=i(k)),\i/j m= 2 (k))} form degenerate pairs at each 
momentum, thus u m =o = ^ 777=3 and Vm—i = v m = 2 - Since 
S 777 = 0, 

tq = ^777=0 = ^777=3 = ^777=1 = ^777=25 ( 20 ) 

hence there is only one independent topological num¬ 
ber, V\. Since the topological charge of a Dirac point 
can be nonzero only when the valence and conduction 
bands have different rotation eigenvalues, a 4-band model 
with Dirac points can be constructed by using a ba¬ 
sis {|^m=o( fc ))> lV4= 3 ( fc )>: l^m= 2 ( fc )>} where 

A = v (c) and B = c (v). Since n^f 0 = n m=z = 
~ n m =1 = ~ n m =2 = ±5 and = —n® for 4-band 
models, the Dirac point has a nonzero topological invari¬ 
ant 

17! = ± 1 . ( 21 ) 


4- Ce symmetric systems 

In the presence of a C 6 rotation symmetry, 
{|^777=0(fe)), |^777=5(&)}}, {|V’m=i(fc)),|V’T 77 = 4 (fc))}, an d 


{\ 4 ’m= 2 (k)), |^777=3(fc))} form degenerate pairs. Thus 
^ 777=0 = ^ 777=1 : ^ 777=45 ^ 777=2 = ^ 777 = 3 - Consid¬ 

ering Vm = dj we can find two independent topolog¬ 
ical numbers (Fi,F 2 ) £ Z , which, for instance, can be 
defined as, 

^1 = ^777—0 = ^777=5: 

1^2 = ^777 = 1 = £'777=4' ( 22 ) 

However, for convenience, we can also use 
(^777=0, ^777=1, ^777=2) to indicate the topological charge 
in which v m —o + Vm=i + Vm =2 = 0. A 4-band model 
can be constructed by choosing two different pairs of 
eigenstates such as 

{|l^777=o)l |^777— 5 ) |V , 777=l)) |f/ , 777=4)}) 

0)7 1^777=5) 1^777 =2); •) (23) 

{|t^777—1)7 1^777=4) 1^777—2)7 1^777—3)}' 

For a given 4-band model, a nonzero topological number 
v m = ±1 can be assigned if J m is the eigenvalue of one of 
the four bands. Whereas v m = 0 if J m is the eigenvalue of 
the other two states which are not included in the 4-band 
model. Therefore the topological charges of the system 
are in the form of 

(^777=07 ^ 777=1 7 ^ 777 = 2 ) = C—C 1 7 “F f- 7 CO 7 

(^777=07 ^777=1 7 ^777=2) = C~^ 1 7 0, = p 1) , (24) 

(^777—0 7 £'777 = 1 , ^ 7U = 2 ^ = (0, A 17 “F 1) 7 

for each case shown in Eq. (PZ51) . respectively. Then the 
corresponding (Fi,F 2 ) are 

(fi,F 2) = (±1, =Fl), 

(F!,F 2 ) = (±1,0), (25) 

(Fi,F 2 ) = (0, ±1), 

respectively. 


Cn 

4-band model 

Materials 

c 2 

Not allowed 


c 3 

Fi = ±1 

NaaBi [2], strained TIN [39j 

c 4 

Fi = ±1 

Cd 3 As 2 [2] 

Ce 

(Fi,F 2 ) = (±1, =Fl) 


C 6 

(Fi,F 2 ) = (±1,0) 


Ce 

(Fi,F 2 ) = (0, ±1) 



TABLE II. Topological charges of class I Dirac semimetals for 
4-band models and relevant materials. 


D. Fermion number doubling in class I Dirac 
semimetals 

Up to now, we have described how to determine the 
topological charge of a single Dirac point. Now let us 
compare the topological charges of two Dirac points at 
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the momenta fco and — fco- Due to the inversion symme¬ 
try satisfying [P, Cn\ = 0, the eigenstates at fc and — k 
satisfy the following relationship, 

C N P\%l)m{k)) = PC N \lj)m{k)) 

= PJm\^m{k)) 

= J m P\^ m (k)), (26) 

which means that if there is an eigenstate with the eigen¬ 
value J m at k , there should be a degenerate eigenstate 
with the same eigenvalue J m at k. This imposes the 
following condition of 

r^’ s (fc 0 ) = n s m N (-k 0 ). (27) 

It is to be noted that the north (south) pole at fco and 
the south (north) pole at —fco are interchanged under the 
inversion symmetry. Thus we obtain 

^m(ko) — l^m( fco). (28) 

Since the net topological charge of the two Dirac points 
related by the inversion symmetry is zero, we obtain the 
following conclusions. 

• The total topological charge of two Dirac points 
within the first Brillouin zone should be zero in 
each angular momentum channel ( J m ), he., 

E^ D) = 0 > ( 29 ) 

I'D 

where id labels Dirac points. It is worth to note 
that this is nothing but the Nielsen-Ninomiya the- 
ore m 36 ' 37 for three-dimensional Dirac semimetals. 

• A stable Dirac point with a nontrivial topological 
charge cannot exist at a time-reversal invariant mo¬ 
mentum (TRIM) where fco = — fco modulo a recip¬ 
rocal lattice vector due to the relationship 

zy m(fco) = ^m( fco) — ^'m(fco) = 0- (30) 

IV. SCREW ROTATIONS AND A SINGLE 
DIRAC POINT 

A. Projective symmetry and circumventing 
fermion number doubling 

It is worth to note that the doubling of the number 
of Dirac points in class I Dirac semimetals results from 
the commutation relation [P, Cn} = 0 as discussed in 
Sec. 1III Dl This means that it may be possible to avoid 
the doubling of the Dirac points, once the commutation 
relation is violated, i.e., [P, Cn] ^ 0. However, the pres¬ 
ence of a single Dirac point on the rotation axis brings 
about a more fundamental problem, when the periodic 
structure of the system is considered. This is because 
the band crossing (or a nonzero topological charge of a 




FIG. 3. Band structure along the rotation axis (z axis) of (a) 
a class I Dirac semimetal and (b) a class II Dirac semimetal. 
J 1,2 and J 3,4 are the rotation eigenvalues of each doubly de¬ 
generate band. A band crossing requires Ji,2 ^ Jza- hr class 
I (II) Dirac semimetals, the band crossing condition and the 
periodicity of the eigenstates are compatible (incompatible) 
when Ji,2,3,4 are constant on the rotation axis. 


Dirac point) requires that the valence band and the con¬ 
duction band should have distinct eigenvalues, whereas 
the lattice periodicity requires the continuity of the eigen¬ 
state and its relevant eigenvalues as described in Fig. [3] 
Therefore the presence of a single Dirac point or an odd 
number of Dirac points on the rotation axis sounds un¬ 
physical, when the rotation symmetry exists along a line 
satisfying the periodic boundary condition. 

One possible way to circumvent the contradiction is 
when the rotation symmetry is realized projectively. 
Namely, if the rotation eigenvalue is well-defined only 
up to an additional phase factor, it is possible to create a 
single Dirac point compatible with the lattice periodicity 
by adjusting the phase degrees of freedom on the rota¬ 
tion axis. In fact, a screw rotation is such an example of 
projective symmetry, which can support a single isolated 
Dirac point as discussed in detail below. 

A screw rotation (Cn, p ) is a non-symmorphic symme¬ 
try operation composed of an ordinary rotation (Cjv) 
followed by a partial lattice translation t p = (p = 
1,..., N — 1) parallel to the rotation axis. Here z is the 
unit lattice translation along the z axis assuming that 
the screw axis is parallel to it. Schematic figures de¬ 
scribing all possible screw rotations in 3D crystals are 
shown in Fig. [4] Let us note that, in many crystals, the 
screw rotation axis does not pass the reference point of 
the point group symmetry, which is invariant under point 
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group operations of the lattice. In this case, the partial 
translation r p associated with the screw rotation CV, P 
also includes in-plane translation components perpendic¬ 
ular to the screw axis direction. Generally, Cn, p can be 
compactly represented as 

C N , P = {C n \t p }, (31) 

where 

T P = ( T p,2u T p,yi T p,z = "77 )■ (32) 


In the real space, CV iP transforms the spatial coordinates 
in the following way, 

C N , P : ( x , y, z) (x' + r PtX , y' + t p ^ v , z + r p>z ), (33) 


where 


2 n 

' TV 

27r 


x = x cos — — y sin —. 

TV ^ AT’ 


y = x sin — + y cos ■ 


27T 

77 : 

27T 


(34) 


iv ' " ■” IV 

Now we consider the combination of a screw rotation 
and the inversion symmetry. At first, we see 

PCn, p • {x, y•) z*) ^ ( x t~p,x > y A 3 ,?/> ^ A *, 2 ) 3 

thus 


PCn,p = {PC N | - r p }. (35) 

Similarly, 

G/V, p P • (*G P) t T A 3 , 2 ’ y ~b "b 3 , 23 ^ d“ A 3 , 2 )’ 


thus 


Cjv iP P = {PCW|t p }. (36) 

Let us note that [P, Cjy] = 0 in general. Equations (1551) 
and (1561) clearly show that generally [P, Cn, p \ ^ 0 due to 
the partial translation t p , thus there is a chance to avoid 
the doubling of the Dirac points. 

In the presence of the screw rotation symmetry C/v, p , 
the Bloch Hamiltonian H(k z ) on the rotation axis (k x = 
k y = 0) satisfies 


C N , p {k z )H(k z )C^ p {k z ) = H(k z ), (37) 

where 

[C N , P (k z )] N = - exp {-ipk z ). (38) 

Here the minus sign stems from the spin 1/2 nature of 
electrons. Therefore all bands on the k z axis can be la¬ 
beled by the eigenvalues of CN, P {k z ) given by 


J m {k z ) = exp 


(2m+1) 

in- 


N 


exp H^ z ) 


= dm exp ( - *|^ 2 ) 3 


(39) 




FIG. 4. Schematic figure describing all possible screw ro¬ 
tations Cn, p in 3D crystals. In each figure, the system is 
periodic under the lattice translation t along the vertical di¬ 
rection. Cn, p indicates a ^ counterclockwise rotation com¬ 
bined with a partial translation j^t. (a) (72,1 symmetry, (b) 
p 3 , P symmetry (p = 1,2). (c) Ca, p symmetry (p = 1,2,3). 
(d) (76 , P symmetry (p = 1, 2, 3,4, 5). 


where J m is an eigenvalue of Cat defined in Eq. (0). It is 
worth to note that the eigenvalue of the screw rotation 
Cn, p is not J m but J m exp(— ijjk z ) which varies along 
the rotation axis. Therefore through the variation of this 
additional phase factor, it may be possible to satisfy the 
condition for the band crossing to create a Dirac point 
and the periodicity (or the continuity) of the eigenval¬ 
ues, simultaneously, even in the presence of a single Dirac 
point. In fact, the assignment of non-quantized quantum 
numbers to fermions, such as J m (k z ) varying in the mo¬ 
mentum space, is one way to get around the fermion dou¬ 
bling problem, as pointed out by Nielsen and Ninomiya 
in their seminal wor k 36 i 37 . 


B. Screw rotations and band connections at the 
zone boundary 


The momentum dependence of screw rotation eigen¬ 
values J m (k z ) shown in Eq. (1391) induces nontrivial band 
connections between different eigenstates at the Brillouin 
zone boundary. For instance, if the system has 27 t peri¬ 
odicity along the k z axis, we find that 


Jm(k z + 2tt) = exp 


(2m 

in - 


— Jm—p{kz)') 


— 2p+l) l 
N 



‘N *') 

(40) 


thus the eigenstate with the eigenvalue J m (k z ) should 
be smoothly connected to the other eigenstate with the 































FIG. 5. The band connection required by a screw rotation 
Cn, p when the system is 2n periodic along the k z direction. 


eigenvalue J m -p{k z ) ( J m +p{k z )) at the Brillouin zone 
boundary k z = tt ( k z = — -jt) as shown in Fig. [5] 
This naturally gives rise to a band crossing point at 
the Brillouin zone boundary. If this band connection 
is compatible with the T and P symmetries, a single 
Dirac point can be realized at the Brillouin zone bound¬ 
ary. The PT symmetry requires that the state with 
the eigenvalue J m (k z ) = J m exp(—ijjk z ) should be lo¬ 
cally degenerate with the other state with the eigen¬ 
value JN-m-i(k z ) = ex p(— i’^kz) at each k z . Sim¬ 
ilarly, we can expect the degeneracy between two states 
with the eigenvalues J m - P (k z ) = J m -p exp(— ijjk z ) and 

JN-m+ P -i(k z ) = Jm- p exp(-i-H-fc z ), respectively. Here 
the important point is that the screw rotation requires 
a nontrivial band connection between Jpf-m-i{k z ) and 
Jiy-m+p-i(kz), similar to the relation shown in Eq. (HUl) . 
Namely, 

JN—m—l{k Z T 27t) = Jn~ m— 1— p{kz) 

= JN—m+p—l[kz)i ( 41 ) 

which gives 

N — m— 1— p = N — m + p — 1 (mod N), (42) 

thus 

N 

P=y (43) 

Since p is an integer, this condition can be satisfied only 
in systems with (72,1, (74,2, C 6,3 symmetries. 

Let us note that, in 3D crystals, the periodicity along 
the k z direction can be longer than 2 ir/a z although the 
system is periodic under the translation by a z along the 
z direction, unless a z is a primitive lattice vector. (For 
instance, it happens in the face centered cubic lattice.) 
Generally, when the system is 2nir periodic along the k z 
axis with an integer 1 < n < N, 


J m {k z + 2mr) = exp 


(2 TO 

ITT - 


— Jm—np(kz)i 


2 np + 1 )' 
N . 



(44) 


thus the eigenstate with the eigenvalue J m (k z ) should 
be smoothly connected to the other eigenstate with the 


eigenvalue J m -np(k z ) ( J m + np (k z )) at the Brillouin zone 
boundary k z = rnr (k z = —mr). Considering the PT 
symmetry and following the same procedure that we have 
used to derive Eq. fl3l) . we obtain 


2 np = 0 (mod N). (45) 


For example, when the system is 47 r periodic (n = 2), 
Eq. (O can also be satisfied in systems with 64 , 1 , Ca, 2 , 
( 74 , 3 ; ( 76,3 symmetries. However, in systems with (7 4,2 
and ( 76,3 symmetries, a 47 r shift merely maps an eigen¬ 
state into itself, hence nontrivial band connection at the 
Brillouin zone boundary is not expected. On the other 
hand, when the system is 67 t periodic (n = 3), Eq. (H5l) 
can be satisfied in systems with (76 , p symmetry where 
p = 1,2, 3,4, 5. However, in the case of ( 76,2 and Ce ,4 
symmetries, a 67t shift connects an eigenstate with it¬ 
self. Also in the ( 76,3 symmetric case, a 67t shift is simply 
equivalent to a 27 t shift, which is already considered be¬ 
fore. Hence only the systems with (76 ,1 and ( 7 e ,5 can sup¬ 
port a nontrivial band connection at the Brillouin zone 
boundary k z = ±37t. 

To sum up, in a Cn, p symmetric system satisfy¬ 
ing p/N = p'/N' with two co-prime numbers p' (an 
odd integer) and N' (an even integer), two distinct 
Cn, p eigenstates should be connected to each other at 
the Brillouin zone boundary k z = ±N'tt/2. Namely, 
the eigenstate with the eigenvalue J m (fc 2 ) should be 
smoothly connected to the other eigenstate with the 
eigenvalue Jm-N/ 2 (k z ) at the Brillouin zone boundary 
k z = ±(iV' , 7 r )/2 in the following way, 


Jm(k z + N'tc) = exp 


(2m — Np' + 1) 
in — 


N 


exp 


Jm—Np'/2(kz) 

Jm—N/2{k , z')i 


( * 

(46) 


where we have used the fact that jf is an odd integer 
and m is well-defined modulo N. It is interesting to note 
that the eigenvalue J m (fc 2 ) at the zone boundary k z = 
±N f 7t/ 2 becomes 


J m {k z = ±N'ir/ 2) = J m exp ( =F *^- 7r ) 

= T*Jm(-l) (p '- 1)/2 . (47) 


Namely, the eigenvalue J m (/c 2 ) is simply given by ±*J m 
at the zone boundary. This additional factor ±i gives 
rise to the following relations between the screw rotation 
Cm n and the inversion P at the zone boundary k+ = 
(0, 6 , ±N'n/2), 


PCn, p 

Cn, p P 


k±) =T*(-l) (p '- 1)/2 PCv|fc±), 

fe±) =±i(-l) (p '- 1)/2 P<7 JV |fc±), 


(48) 
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FIG. 6 . (a) An example of the band connection required by 

the screw rotation in Eq. m- For each band, the eigenvalue 
of the screw rotation 62,1 is marked in the figure. (b,c) Pos¬ 
sible band structure protected by a two-fold screw rotation 
C 24 with a Dirac point at k z = 0 and at k z = n, respectively. 

which can be easily derived from Eq. (1551) and (1351) . 
Hence P and Cn, p anticommute when the Bloch state 
| k± = (0,0, ±N'n/2)) is used as a basis for the repre¬ 
sentation. From this, we obtain the following general 
principle to create a stable Dirac semimetal with a single 
Dirac point. Namely, 

• The Dirac point should be located at a TRIM at 
the Brillouin zone boundary ( k z = ±N'tt/2) where 
the screw rotation symmetry anti-commutes with 
the inversion symmetry. 


C. Applications 

In the following, we examine the possible Dirac 
semimetals with a single Dirac point by considering var¬ 
ious screw rotation symmetries explicitly. 



FIG. 7. An example of the band connection required by 63,1 
symmetry shown in Eg. (1541) . Two bands with the eigenvalues 
Jo(k z ) and Ji(k z ) ( Jo(k z ) and Ji(k z )) form a conduction (va¬ 
lence) band. However, each degenerate pair violate the PT 
symmetry. 


Now we prepare two bands U/o (k z ) and dti (k z ) with an 
eigenvalue Jo(k z ) and Ji(fc 2 ), respectively, and construct 
a band structure with a Dirac point. Here the cru¬ 
cial point is that the band d>o(A: 2 ) (Ti (k z )) should be 
smoothly connected to the other band ’Fi (k z ) (’Fo(fcz)) at 
the Brillouin zone boundary ( k z = ± 7 r) to satisfy Eq. (1501) 
and (I5T1) as shown in Fig. [G] This naturally gives rise to 
a band structure with a single band crossing point at a 
TRIM. Considering the P or T symmetry, there are two 
possible band structures having a single band crossing 
point as shown in Fig. [5] (b) and (c). In each case, the 
band crossing point locates at a TRIM either at k z = 0 
(Fig. [G] (b)) or at k z = tt (Fig.[5](c)). However, let us note 
that, due to the TP symmetry, the state with the eigen¬ 
value Jo{k z ) (Ji(k z )) should be locally degenerate with 
the other state with the eigenvalue J\{k z ) ( Jo(k z )) at 
each momentum k z . This requires that both conduction 
band and valence band should have the same eigenvalues 
of Jo(k z ) and Ji(k z ), hence a stable Dirac point cannot 
be created due to the finite hybridization between the 
valence and conduction bands. 


1. Two-fold screw rotation C 2 ,i 


2. Three-fold screw rotation 


A two-fold screw rotation symmetry ( 72,1 has the fol¬ 
lowing two eigenvalues 


Jo 

Ji 


( k z ) = J 0 exp ( - *2 kz ) = exp 
(k z ) = Ji exp ( - = exp 


% 2 {kz 
*2 {kz 


oj. 

3?r) , 
(49) 


where 


Jm(k z — 7f) ^ Jm^kz; — 7r), (h0) 

for m = 0,1 and 

Joikz) = Ji(fcz + 2 tt), 

Ji(k z ) = J 0 (k z + 2n). (51) 


In the case of a three-fold screw rotation < 73 , 1 , there 
are three possible eigenvalues given by 


Jo 

Ji 

Ji 


( k z ) = J 0 exp - *gfczj = exp 
( k z ) = Ji exp ( - *^z) = exp 
(k z ) = Ji exp ( - *^z) = exp 


i 3 (k - 

'b k ‘ 

'b k ‘ 


3ti-) , 

5tt) , 
(52) 


where 


Jm{k'Z — 7r) ^ — Tr)? 


(53) 
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for m = 0 , 1,2 and 

MM = Ji{k z + 27r), 

MM = J 2 (k z + 2rr), (54) 

MM = J 0 (fc 2 + 27r), 

Now we prepare three bands 4/0,1,2(M with an eigen¬ 
value Jo, 1,2 (fez) > respectively, and construct a band struc¬ 
ture with a Dirac point. To satisfy Eq. (l53ll and (l54ll . each 
band with a given 63,1 eigenvalue should be connected 
to the other two bands with different 63,1 eigenvalues at 
each Brillouin zone boundary. An example of the band 
connection satisfying Eq. (15H) is shown in Fig. [7] Ac¬ 
cording to Fig. [3 when the band structure is drawn for a 
reduced Brillouin zone with k z £ [—7r, 7 t], two bands with 
the eigenvalues Jo(M and Ji(k z ) (or J 0 (M and J 2 (fc 2 )) 
would form a degenerate band. 

However, this band structure is generally incompatible 
with the T and P symmetries. The PT symmetry re¬ 
quires the band with an eigenvalue Jo(fc 2 ) to be locally 
degenerate with the other band with an eigenvalue J 2 (fc 2 ) 
whereas the band with an eigenvalue J\(k z ) to be locally 
degenerate with the other band with the same eigenvalue 
J\{k z ). Namely, the four bands 

{M {k z ), MM; MM, MM} ( 55 ) 

would form a basis to create a Dirac point. In the paren¬ 
thesis, the first two bands form a conduction (valence) 
band whereas the last two bands form a valence (con¬ 
duction) band. This basis is obviously incompatible with 
the band connection described in Fig. 0 It is straight¬ 
forward to show that the same problem happens for C 32 
symmetric systems. Therefore a system with a three¬ 
fold screw rotation cannot satisfy the PT symmetry at 
the same time, hence cannot have a stable Dirac point. 
In fact, every screw symmetric system which does not 
satisfy Eq. (1751) or Eq. (H5l) has the same problem, thus 
a Dirac semimetal with a single Dirac point cannot be 
created. 


^ 0 , 3 ^ z ) 



FIG. 8 . An example of the band connection required by C. 4,1 
symmetry shown in Eg. (1601) . It is assumed that the system 
is 47 t periodic. The band with the eigenvalue Jo{k z ) (Js(k z )) 
should be connected with the other hand with the eigenvalue 
J 2 (fc 2 ) (Ji(fc 2 )) at the zone boundary. Moreover, due to the 
PT symmetry, two bands with the eigenvalues Jo(k z ) and 
J 3 {k z ) (or, J 2 (fc z ) and Ji(k z )) should be locally degenerate at 
each momentum. In this case, the band connection required 
by the screw rotation is compatible with the PT symmetry. 


If the system is 27 t periodic along the k z direction, we 
obtain 

Jm(k Z — Tf) 7 ^ Jm[kz — 7r); (57) 

for m = 0,1, 2, 3 and 

Jm(k z ) = Jm+i(k z + 27t), (58) 

hence the state with the eigenvalue J m (M should be con¬ 
nected to the state with the eigenvalue J m +i(fc 2 ) at the 
Brillouin zone boundary ( k z = ± 7 r). However, this band 
connection is not compatible with the PT symmetry of 
the system, which requires the state with Jo(fc 2 ) ( Ji(k z )) 
to be degenerate with the state with Js(fc 2 ) (J 2 (fc 2 )). 

On the other hand, if the system is 47 t periodic along 
the k z direction, 


3. Four-fold screw rotation C 4 ,i 


Jm=o, 1 , 2 , 3 (k z — 27t) zfz J m= o, 1 , 2 , s(k z — 27t), (59) 


In the case of the four-fold screw rotation C 4 ,i, there and 
are four possible eigenvalues given by 


Jo(M = Jo exp ( 

M) 

= exp 

Jl(M = Ji ex P ( 


= exp 

MM = ^2 exp ( 


= exp 

MM = ^3 exp ( 

,M) 

= exp 


i\{k z 

l ~(k z 

i\(k z 

i\(kz 


t0|. 

3tt) , 

5tt) , 

7ir) . 
(56) 


Jm{k z ) = j m+2 (fc 2 + 47 t). (60) 


This band connection is compatible with the PT symme¬ 
try. The basis for creation of a Dirac semimetal with a 
single Dirac point at the Brillouin zone boundary is given 
by 


{*o(fc 2 ), 'MM; MM, MM}. 


(61) 
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4- Four-fold screw rotation C.4,2 


6. Six-fold screw rotation Ce ,3 


In the case of the four-fold screw rotation C 4 2 , there 
are four possible eigenvalues given by 


Jo{k z ) = Jo exp ( - k z ) = exp 


~ / 1 \ 

1, 3 1 

J 1 (k z ) = Ji exp - i-k z J = exp 

[-*2^-2^] 

~ / 1 \ 

1 5 j 

J 2 {k z ) = J 2 exp - *-fc 2 J = exp 

[~ l -(kz~-n)\ 

~ / 1 \ 

1 7 i 

J 3 \k z ) = J 3 exp - i 2 kz ) = exp 



(62) 

where 

Jm{k z = - 7 r) ± Jm(k z = 7 r), (63) 

for m = 0 , 1 , 2 ,3 and 

Jm{kz ) = Jm+2{kz T 27t), 

hence the state with the eigenvalue Jo(fc 2 ) (Ji(fc z )) 
should be connected with the state with the eigenvalue 
J 2 {kz ) {J:i{k z )) at the Brillouin zone boundary (k z = 
± 7 r). This band connection is compatible with the PT 
symmetry. The basis for creation of a Dirac semimetal 
with a single Dirac point at the Brillouin zone boundary 
is given by 

{*o (k z ),* 3 (kz); *i(fc„), ^ 2 (k z )}. (64) 


In the case of the six-fold screw rotation Ce t 3 , there 
are six possible eigenvalues given by 


J r n (/c 2 ) — J m exp ^ i ^ k z ^ 


= exp 


2 m + 1 


(69) 


where m = 1,2, 6 . When the system is 2 - 7 T periodic 
along the k z direction, these eigenvalues satisfy 


Jmik, = 7r) ^ Jm(k z = 7r), (70) 

and 

Jm(k z ) = Jm+z{k z + 27r), (71) 


hence the state with the eigenvalue Jo,i, 2 (k z ) should be 
connected to the state with the eigenvalue < 73 , 4,5 (fc 2 ) at 
the Brillouin zone boundary ( k z = ±7r), respectively. 
Due to the T and the P symmetries, the pair of states 
with the ^eigenvalues {J 0 (fc 2 ): ^(M}, {Ji(k z ), Ji{k z )}, 
{<7 2 (fc 2 ), J 3 (k z )} should be locally degenerate at each mo¬ 
mentum k z . Considering the band connection described 
in Eq. (EH we find the following basis 

{*o(M, ^5 (kz)-, * 3 (*,), tf 2 (fcz)}, (72) 

which can create a Dirac semimetal with a single Dirac 
point at the Brillouin zone boundary. 


5. Four-fold screw rotation 64,3 


7. Six-fold screw rotation Ce :P ^3 


In the case of the four-fold screw rotation 64 , 3 , there 
are four possible eigenvalues given by 


Jo{k z ) = Jo exp ( - = exp 

~i\{kz - ^tt) 

~ / 3 \ 

3, 3 j 

J 1 (k z ) = Ji exp - i-fc 2 J = exp 

-i-^kz - - 7 r) 

~ / 3 \ 

3, 5 y 

J 2 {k z ) = J 2 exp - i-fc 2 J = exp 

~ l ^k z - g7r) 

~ / 3 \ 

3, 7 ,1 

7 3 (Kz) = J 3 exp - i-k z J = exp 

k z - g7r) 


Similar to the case of C& t \ symmetric systems, a Dirac 
semimetal with a single Dirac point can be created only 
if the system is 47 r periodic along the k z direction. Then 

Jm{k z = -27t) zfz J m (k z = 27t), (66) 


In the case of the six-fold screw rotation Cq jP ^ 3 , a sin¬ 
gle Dirac point can be created only if the system is 67 t 
periodic along the k z direction. 

In ( 76,1 symmetric systems, there are six possible eigen¬ 
values given by 

Jm(kz) = Jm exp ^ i~k z ^j 

= exp | — i^[k z — (2m + l)7r]|, (73) 

where m = 1,2,..., 6 . These eigenvalues satisfy 

Jm {kz = -37r) ^ J m {k z = 37r), (74) 

and 

Jm{k z ) = Jm+3 ik z + 6tt), (75) 


for m = 0 , 1 , 2 ,3 and 

Jm{k z ) = Jm- 2 {k z + 47t). (67) 

The basis for creation of a Dirac semimetal with a single 
Dirac point at the Brillouin zone boundary is again 

{^o(fcz), ’MM; ^1 {k z ), 'f> 2 {k z )}. (68) 


In ( 76,5 symmetric systems there are six possible eigen¬ 
values given by 


Jm[k Z ) — Jm exp ^ \kz^ 

= exp 


r i 5 ( k 

2 TO +1 \ ' 

[ l 6\ kz 

5 ' /J 


(76) 
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where to = 1,2, 6 . These eigenvalues satisfy 

J m (k z =-3n) ^ J m (k z = 3n), (77) 

and 

Jm{k z ) = Jm+3(kz + 67t). (78) 

In both C 6 ,i and ( 76,5 symmetric cases, we find the fol¬ 
lowing basis 

{*o(fc*)> ^5 (k z ); $3 (*,), ^ 2 (fc 2 )}, (79) 

which can create a Dirac semimetal with a single Dirac 
point at the Brillouin zone boundary. 

Finally, it is straightforward to show that ( 76,2 and 
( 76,4 symmetric cases, which are similar to the system 
with a three-fold screw symmetry, cannot support a Dirac 
semimetal with a single Dirac point. 

V. CLASS II DIRAC SEMIMETALS 

Based on the discussion in Sec. IIVI we define class II 
Dirac semimetals in the following way. Class II Dirac 
semimetals are associated with a special type of rotation 
symmetry Cn which anti-commutes with the inversion 
symmetry, i.e., 

{P,C N } = 0. (80) 

Considering that i) a rotation operator (7 generally has 
a form of C = exp(*@J) with the angular momentum 
operator J and the rotation angle 9 , and ii) J is a pseu¬ 
dovector which is even under the inversion symmetry P, 
i.e., PJP -1 = J, the anti-commutation relation between 
P and Cn looks quite unusual. However, as discussed 
in the previous section, non-symmorphic screw rotation 
symmetries can generally satisfy such anti-commutation 
relations at the Brillouin zone boundar y 27 ' 28 , when Bloch 
states are used as a basis for the representation of P and 
Cn- We first describe the physical consequence resulting 
from the relation {P, Cn} — 0 and the topological charge 
of the associated Dirac point. After that we describe how 
the Dirac point protected by a screw rotation symmetry 
leads to a class II Dirac semimetal, and its associated 
topological charge can be determined. 

A. Symmetry constraint and band connection 

Let us consider the constraints to the rotation eigen¬ 
values J m due to the symmetries P, T, Cn satisfying 

[T,C n ]= 0 , {P 1 C N } = 0, [T,P]=0, (81) 

(82) 


Moreover, considering Eq. dUD, we assume the Cn eigen¬ 
value J m to have the following form, 

Jm — iJmi (83) 

where J m = exp[i 7 r( 2 m + 1 )/N}. The action of the PT 
on an eigenvector \i/j m (k)) of Cn with an eigenvalue J m 
gives 

C N PT\^ m (k)) = -PTC N \^ m (k)) 

= -PTJ m \il> m (k)) 

= -r m PT\^ m {k)) 

= J N -m-lPT\lp m (k)) 

= Jn~ m — 1 j IpN—m — 1 (A?)) 5 (84) 

from which we find PT\ip m ) is a degenerate eigenvector 
with the eigenvalue JN-m-i- On the other hand, P and 
T transform the eigenvector \tp m ) at k to another eigen¬ 
vector at — k as 

C N T\^ m (k)) = TC N \il> m {k)) 

= T Jm\l/j m (k)) 

= J* m T\^m(k)) 

= J^ m -lT\^m{k)) 

= JN_ m _ 1 \ipN_ m _ 1 (-k)), (85) 

and 

C N P\ipm(k)) = -PC N \i>m{k )) 

= -PJm\4>m(k)) 

= -J m P\ip m {k)) 

= Jj* +m P\lpm(k)) 

= JN +m \ipN +m (-k)). (86) 

Since y — to — 1 and -^ + to should be an integer, we 
find that N should be an even number. Hence the class 
II Dirac semimetal cannot exist in systems with C 3 sym¬ 
metry. From Eqs. (|84ll . ([85jl. and (l 86 l) . we find that if 
\ip m (k)) is an eigenstate at fc, \%})N-m-i{k)) is also a de¬ 
generate eigenstate at the same momentum k, whereas 
k)) and k)) are degenerate at — k 

with the same energy as \ip m (k)). From this, we can infer 
the band connection near a TRIM (/strim) where fc and 
k are equivalent. The question is whether two doubly- 
degenerate states at k and — k are crossing or smoothly 
connected at Antrim • A smooth connection of degenerate 
bands requires the rotation eigenvalues of the states at k 
and — fc to be identical, which is obviously not satisfied 
in this case. Therefore there should be a band crossing 
at A;tr,im where the Dirac point locates. 


T 2 = -1, P 2 = 1. 
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(a) N = 4n (b) N = 4n+2 



FIG. 9. Constraints on the rotation eigenvalues (a) for TV = 
4n, (b) for TV = 4n + 2 with an integer n. The (black) solid 
circle indicates a unit circle in the complex plane, and each 
(red) dot on the circle denotes J m . Four dots connected by 
a dotted line are related by the P, T, and PT symmetries, 
hence only one of them is independent. In the case of (b), the 
topological charge associated with J(jv-2)/4 and J( 3 n-2)/4 is 
zero. 

B. Topological charge 

From Eqs. I(M|) . (1551) . and (1551) . we can easily find 
the constraints to the topological charge v m of the Dirac 
point locating at a TRIM. At first, the PT symmetry 
requires that 

Vf-m-l = U f+m- ( 87 ) 

Moreover, since P or T symmetry interchanges the 
north and south pole surrounding a Dirac point at a 
TRIM, we find 

V m = VM-m-l = _ m _i = ~VN +m . (88) 

This constraint reduces the number of independent topo¬ 
logical numbers i/ m . Considering that TV is an even in¬ 
teger, we distinguish two cases, i.e., when TV = 4n and 
when TV = 4n+2 with an integer n. In each case, the rela¬ 
tion between different eigenstates are described in Fig. [51 
from which we find the topological invariant 

N 

(r-o,..., ^n_ 1 ) g if TV = An, 

(vn .z/ jv -6 ) g Z 4 if TV = 4n + 2. (89) 

4 

Hence the topological charge of the system with C 4 or Cq 
symmetry is an element of Z whereas C 2 or C 3 symmet¬ 
ric system cannot support a Dirac point with a nonzero 
topological charge. 


Cn 

Topological charge 

c 2 

Not allowed 

c 3 

Not allowed 

Ci 

Z 

Ce 

Z 


TABLE III. Summary of topological charges of class II Dirac 
semimetals. 

pairs at each momentum fc, and these four states cross 
at a TRIM, and create a Dirac point. A 4-band model 
can be constructed by using these four states. 


1. C 2 symmetric systems 

Possible J m values are J m - 0 = *exp(*i7r) and 

J m= 1 = iexp(*|7r). Due to the PT symmetry, 

{|V’m=o(fe)),|V , m=l(fc)}} and {|V’m=o( fe ))>IV’m=l( fe )>} 
form degenerate pairs. The symmetry constraint in 
Eq. (1551) requires 

l / m =0 = 1 = 0 . ( 90 ) 

Hence a C 2 invariant system cannot support a stable 
Dirac point at a TRIAL 


2. C 3 symmetric systems 

Possible J m values are J m =o = iexp(*i7r), J m= i = 
iexp(f-7r) = —i, and J m =2 = *exp(i|7r). In the case 
of |^ m =i), the P symmetry requires that CzP\il) m -\) = 
-PC 3 |^ m= i) = iP\if> m =i). Thus P\ipm=i) should be an 
eigenstate of C 3 with the eigenvalue +i, which is not 
allowed. Hence a C 3 invariant system cannot support a 
stable Dirac point at a TRIM. 

3. Ca symmetric systems 

Due to the PT symmetry, {|V’™=o), |^m= 3 )} and 
{|^m=i)> |V , m= 2 )} form degenerate pairs. The constraint 
in Eq. (155)) requires that 

^1 — Vm =0 = ^ 771—3 = ^ 771=1 = ^ 777 — 2 ) ( 91 ) 

thus there is only one independent topological invariant 

Vi = ± 1 . 


C. Applications: classification of stable Dirac 

points in 4-band systems 4- Ce symmetric systems 

Due to the PT symmetry, \ipN_ m _i(k))} Due to the P and T symmetries, we find a ba- 

and {\it’N/2-m-i(k)),\4> N /2+m{k))} form degenerate sis {|^ m=0 ), |^ m=5 ); 1 ^= 2 ), |'0-m.=3)}- The constraint in 
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^1 — l / m=0 — On=5 — 2 — ^m= 3? (92) 

thus there is only one independent topological invariant 

Vi = ± 1 . 


Cn 

4-band model 

Materials 

Ci 

Not allowed 


c 3 

Not allowed 


c 4 

Vi = ±1 

b-B10 2 [25] 

c 6 

Vi = ±1 



TABLE IV. Topological charges of class II Dirac semimetals 
for 4-band models. 


D. Topological charge of Dirac points protected by 
screw rotations 


Here we show that the Dirac semimetals protected by 
screw rotations belong to the class II, thus the topological 
charge of the relevant Dirac point can be determined by 
following the prescription described in Sec. IV Bl In partic¬ 
ular, we resolve a subtle issue associated with the multi¬ 
valued nature of the screw rotation eigenvalues, which we 
encounter when we define the topological charge of the 
Dirac point at the Brillouin zone boundary. To under¬ 
stand this, let us again introduce a sphere in the momen¬ 
tum space surrounding the Dirac point, and consider the 
two points fcjv and fcg on the sphere passing the rotation 
axis. To compare the zero-dimensional topological num¬ 
bers at these two points, we need a single-valued wave 
function which varies smoothly around the Dirac point. 
However, since the eigenvalue of a screw rotation J m (k z ) 
is multi-valued, the relevant eigenstates also change dis- 
continuously at the Brillouin zone boundary. To remedy 
this problem, we propose a way to construct a smooth 
function which is single-valued around the Dirac point 
by modifying eigenvectors of Cn, p - Moreover, we show 
that such a smooth single-valued function satisfies the al¬ 
gebraic relations shown in Eq. (1511) . (1551) . and (1551) . thus 
we prove that a Dirac semimetal protected by a screw 
rotation belongs to the class II. 

First, we suppose that p/N = p'/N' holds with an even 
integer N' and an odd integer p' that are coprime, and 
the Brillouin zone boundary is located at k z = N'ir/2. To 
construct a smooth single-valued wave function around 
the Brillouin zone boundary at k z = N'n/2, we prepare 
two eigenstates |H , m (fc 2 )) and |^ m _Ni_ p (k z )) with Cn, p 
eigenvalues J m {k z ) and J m _Nj_ p {k z ), respectively. Then 
we define a hybrid wave function around the zone bound- 


\ipm(8k z )) 


|'3 / m (iV , 7r/2 + 5k z )), (Sk z < 0) 

l^ m _ Jv: p (-iV , 7r/2 + Sk z )), ( Sk z > 0) 

(93) 


which is smooth and single-valued around the zone 
boundary at 5k z = 0; see Eq. (1151) . It is straightfor¬ 
ward to show that |«/> m (^)) is also an eigenvector of 
Cn, p satisfying 

CN, P \i>m{5k z )) = J' m {5k z )\i> m (5k z )), (94) 


where 

J L( Sk z) = exp 


2 to +1 p' 

~N ~2 


exp 


V : 


(95) 


Now we determine the Cn, p eigenvalue of T\ip m (5k z )). 
At first, if 6 k z < 0, we obtain 

T$ m (5k z )) = T\* m (N'n/2 + Sk z )) 

oc |'FAr_ m _i(-Al , 7r/2 - Sk z )) 

= $ N _ m _ 1+ %. p (-6k x )), ( 96 ) 

where we have used the following relation, 

r m = J N -m- 1. ( 97 ) 

On the other hand, if Sk z > 0, 

T$m{Sk z )) = T\* m _^ p (-N'Tr/2 + Sk z )) 
oc \* N _ m _ 1+ q. p (N'ir/2 - Sk z )) 

= |^W_ m _ 1+ 3i p (Sk x )). (98) 

From Eqs. (1151) and (1551) . we obtain 

C N , p T$ m (5k z )) = J' N _ m _ 1+ ^ p (-6k z m m (6k z )) 

= J , N - m - 1+ %(Sk ! .)T\ : ijj m (6k z )), 

= J'%- m - 1 {-8k z m m (Sk z )), ( 99 ) 

where we have used 

N' N N 

-yP = - jp ' = — ( mod N ), ( 10 °) 

in which p' is an odd integer. 

The Cn, p eigenvalue of P\i() m (Sk z )) can also be ob¬ 
tained similarly from the following two relations, 

P$m(6k z < 0)) = P\^ m (N'TT/2 + 5k z )) 

(X \'fr m (—N , Tr/2 - 5k z )) 

= 

= | Tp m+ »(Sk z )), (101) 
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and 


P\i/> m (6k z > 0)} = P\* m _^ p (-N'n/2 + 5k~)) 
oc \'& rn _N± p (N'TT/2 - Sk z )) 

= I $ m -%. p (-8k z )) 

= $ m+ »(-6k z )), (102) 

where we have used the fact that m is well-defined mod¬ 
ulo N. From Eqs. (11011) and (11021) . we obtain 

C N , p P$m{Sk z )) = J' m+f {-Sk z )P\i’ m {Sk z )}. (103) 

Finally, the Cn, p eigenvalue of PT\i/j rn (Sk z )) can also 
be obtained by following similar steps, which give 

C N , P PT\4> m {5k z )) = f N - m -i{5k z )PT$ m (5k z )). (104) 

These transformation laws should be compared with 
Eqs. (glD, (1551) . and (1551) . which differ from Eqs. (1M1) . 
(11031) , and (11041) merely due to the momentum-dependent 
phase factor in J^^Skz). However, since the wave func¬ 
tion |i pm(Sk z )} and its associated the eigenvalue J' m (5k z ) 
are smooth and single-valued around Sk z = 0, i.e., near 
the Brillouin zone boundary, the topological charge of 
the Dirac point can be defined by using J' m {5k z = 0) at 
the two momenta k n and k$ near the Dirac point. Once 
5k z is fixed to be 6k z = 0, Eqs. (l84l) . ([85]) . and (l86l) are 
identical to Eqs. (1M1) . (11031) . and (11041) . which shows that 
the Dirac point at the Brillouin zone boundary protected 
by a screw rotation belongs to the class II. 


E. Example 1: a class II Dirac semimetal on a hep 
lattice 

To illustrate the role of screw rotations on the pro¬ 
tection of a Dirac point, let us consider a tight-binding 
Hamiltonian on a hexagonal close-packed (hep) lattice, 
which corresponds to the space group PQ^/mmc (no. 
194). The hep lattice is generated by the primitive lat¬ 
tice vectors ai = ax, a 2 = a(fi + y ), a 3 = cz, 
and two sites in a unit cell located at rj = 0 and 
r 2 = |ai + |a 2 + fa 3 , respectively. The crystal has 
a 6 -fold screw rotation symmetry ( 76,3 about the z axis 
centered at |ai + ^a 2 accompanied by a partial transla¬ 
tion -ici 3 as shown in Fig. 1101 To confirm the presence of 
a single Dirac point at the zone boundary, corresponding 
to the A point in Fig. [11] (b), we construct the following 
tight binding model, 

H = ~ *1 °i C 0 “ *2 C I T * C j 

(ij) (ij) 

+ zA SO V ij c i a zT z Cj + lA SO V ij c i a z T z c Ji 

(ij) ((ij)) 

( 105 ) 





FIG. 10. (a) Hexagonal close-packed (hep) lattice structure 

and 6-fold screw rotation. Two sublattice sites are marked 
with different colors. The arrows indicate a 6-fold screw rota¬ 
tion about the 2 axis (C 6 , 3 ). The number on a lattice site sym¬ 
bol indicates its z-coordinate in the unit of the vertical lattice 
spacing c. (b) Projection of the lattice to the xy plane, (c) 
Spin-dependent complex hopping process between the near¬ 
est neighbor sites (A^q) and the next nearest neighbor sites 
(Agg). When a spin-up electron on the A sublattice hops par¬ 
allel (anti-parallel) to the arrow direction, the corresponding 
Vij = +1 (Vij = - 1 ). 


where t\ (< 2 ) indicates the nearest-neighbor hopping be¬ 
tween the same (different) sublattice sites and AgQ (A^q) 
denotes the spin-orbit interaction between the nearest- 
neighbor (next-nearest-neighbor) sites, 
is +1 (— 1 ) if the bond ij is parallel (anti-parallel) to 
the arrow direction on the bond as shown in Fig. [TO] (c). 
&x,y,z (x x , y , z ) are Pauli matrices indicating the spin (sub¬ 
lattice) degrees of freedom. 

In the momentum space, the Hamiltonian becomes 


H{ k) = F 0 + 


(if } 


so 


p( 2 )j cr 

^so ) az 


Ffl 2 


~(Fi 


FA 

(i) 


1 - 1-2 


so + F^)a z 


where I 2 indicates a 2 x 2 identity matrix and Fop and 
FgQ 2) are given by 

F 0 = —2ti {cos(fc • ai) + cos(fc • a 2 ) + cos[fc • (ai — < 12 )]} , 
Fi = — 2 t 2 cos (e ik bl + e lk b2 + e lk ba ) , 

FgQ = 2AgQ {sin(fc • a±) — sin(fe • a 2 ) — sin[fc • (ai — a 2 )]} , 

FgQ = — 2AgQ [sin(3fe • bi) + sin(3fc • b 2 ) + sin(3fc ■ 63)], 

(106) 


where 61 = §x + ^y, b 2 = -§£ + ^y, and 63 = 
—j%y- We choose t\ = 1, t 2 = 5, and AgQ = AgQ = 5. 
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FIG. 11. (a) Band structure of the tight binding model on a 

hep lattice. There are two Dirac points at the A and L points, 
(b) First Brillouin zone of the hep lattice. 


The resulting band structure is shown in Fig. [TT] We 
can clearly see two Dirac points at A and L , respectively. 
Between these two points, the Dirac point at A is the one 
protected by the 6-fold screw rotation, hence is located 
at the boundary of the rotation axis (z-axis). 

To confirm the symmetry protection of the Dirac point 
at A and its characteristic dispersion, let us examine the 
symmetry of the Hamiltonian. The symmetries, which 
are important for the protection of Dirac points, are the 
time-reversal T, the inversion P, the 6-fold screw rotation 
^ 6,3 = {Ce\§z}, and the glide symmetry M y = {M y ||z} 
where M y transforms the spatial coordinate (x, y , z) to 
(x, — y, z). To find the matrix representation of each sym¬ 
metry operator, we can use the following information. At 
first, for k —> k, we find 


Fo(-fc) = F 0 (k), 

ReFi(-k) = RePi(fc), 

ImFi(— k) = — ImFi(fc), 

^so 2) (-fe) = -^so 2) (fc)> (107) 

which gives 

P = t x , T = ia y K , (108) 


where K is a complex conjugation operator. Moreover, 
under 7 t/ 3 rotation about the z axis, we obtain 


( k x + iky ) -A ( k' x + ik' y ) = ( k x + ik y ) exp 

k z —>■ k z = k z 


and 


F 0 {k') = F 0 (k), 
ReFi(k') = RePi(fc), 
ImPi(fc / ) = — ImPi(fc), 

F^ 0 2 \k')=-F^ 2 \k) 1 


thus 

C e , 3 (k z ) = t x exp exp 


(109) 


( 110 ) 

(111) 


where Ce ) 3 (k z ) means the representation of 65,3 in a 
Bloch basis, in which the momentum dependent phase 
factor results from the partial lattice translation along 
the z direction. Finally, for k y —> —k y , we find 

Po (Me, k y , k z ) — Fq ( k x , k y , k z ), 

RePi(/ca,, k yj k z ) — RePi {k x , k y , k z ) : 

ky,k z ) — luiFi {k x , k y , k z ), 

F so (Mb. ~ k v k *) = F so ( k x,k y , k z ), 

F so (Mr, ~ k v Mr) ~ ~ F so (Mr, ky, kz), ( 112 ) 


thus 


My(k z ) = ir x (jy exp 



(113) 


Let us note that Fg 0 ' term breaks the glide mirror M y . 

On the k z axis with k x = k y = 0, the Hamiltonian 
becomes 


H(k z )=F 0 +F 1 (k z )T x , (114) 


from which we find two degenerate eigenstates, 



(A 


(o\ 

'*»> - 71 

0 

1 


1 

0 


\o) 


W 


with the eigenvalue E + (k z ) = F 0 + F 3 (k z ), and the other 
two degenerate eigenstates, 



( 1 ^ 


( 0 \ 

1 

0 

, , . 1 

1 

71 

-1 

■ w - 2> =75 

0 


1 0 ) 


l-i / 


with the eigenvalue F_(fc z ) = Fq — F\(k z ). Let us note 
that \tp+i), 1 ^+ 2 ), I■0— 1 ), 1 ^- 2 ) are also the eigenstates 
of (76,3 with the corresponding eigenvalues exp(—+ 
if), exp(-i^ + i^-), exp(-i^ + exp(-i^p- + 
i ^), respectively. They are exactly the C 6 ,3 eigenstates 
[See Eq. (1721) ]. which can support a single Dirac point 
at the zone boundary on the rotation axis, i.e., at the A 
point with the momentum k = (0, 0, y-). The low-energy 
Hamiltonian near the A point is given by 


Ha (<?) ~ 3 t 2 q z 1 - 


Qx + 9 y 
12 


T x + 


(2 


24-\/3 


( 3 q 2 x q y - ql)q z r , 


+ 


3^, 


U 


^so( 3 fe 9 y - ql) + —T^io^qlqy - q y ) 


T Z (T Z 


(115) 


where the momentum q are measured relative to the A 
point assuming a = c = 1 and the constant term Fq is 
dropped. It is interesting to note that the dispersion on 
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the ( q x ,q y ) plane is cubic whereas it is linear along the 
q z direction. The cubic dispersion arises due to the an¬ 
gular momentum difference between the conduction and 
valence bands, as indicated in Eq. CZ2D. Thus we obtain 
a cubic Dirac point at the A point, which is protected by 
the 6-fold screw rotation Ce, 3 - 

Finally, let us briefly explain the physical origin of 
the Dirac point at the L point with the momentum 
k = (0, At the L point, the system is invariant 

under a set of point group symmetry operations, which 
is so-called the little co-group at L, G . The little co¬ 
group G is generated by three symmetry operations, a 
glide mirror M y , the inversion P and the two-fold ro¬ 
tation about the y axis C 2 y ■ Due to the partial lattice 
translation t = involved in M y , M y and P do not 
commute. Namely, we find that 

PMy : (x,y,z) (~x,y,-z- , (116) 

thus 

PM y = [PMy I - PJ . (117) 

On the other hand, 



FIG. 12. (a) Band structure of the Fu-Kane-Mele model on 

a diamond lattice. There is a Dirac point at the X point, (b) 
First Brillouin zone of a diamond lattice. 



M y P :(x,y,z)-> (-x,y,-z+^j , (118) 

thus 

M y P = {PM y \Zz} . (119) 

Therefore M y and P anti-commute at the zone boundary 
with k z = 7r/c, which guarantees the four-fold degener¬ 
acy. Since two-fold rotation symmetry cannot support a 
stable Dirac point, C 2y symmetry cannot play an impor- 
tant role here. Moreover, the Fg 0 term breaking the M y 
vanishes at the L point. To establish a general theory 
about the protection of a Dirac point by a glide mirror 
symmetry and its associated topological charge is an in¬ 
teresting research topic, which we leave for future study. 


F. Example 2: a class II Dirac semimetal on a 
diamond lattice 

As a second example of class II Dirac semimetal, let 
us consider the Fu-Kane-Mele Hamiltonian on a diamond 
latticed, 

H = tJ2 4 c i + 8i ^r X c i ff ' ( d lj x d i,) c ^ 120 ) 

(ij) ((«» 

where the first term indicates the nearest-neighbor hop¬ 
ping and the second term connects the second-nearest- 
neighbors with a spin dependent amplitude, rfx are 
the two nearest-neighbor bond vectors traversed between 
sites i and j, and cr X} y }Z are Pauli matrices indicating the 


FIG. 13. Structure of a diamond lattice. Two sublattice sites 
are marked by using different colors. The arrows indicate a 
4-fold screw rotation about the z axis. 


spin degrees of freedom, a denotes the cubic lattice con¬ 
stant. In the momentum space, the Hamiltonian becomes 


H( k) 


£?=i Fi°i *ol 2 \ 

f 0 i 2 


where I 2 indicates a 2 x 2 identity matrix and To, 1 , 2,3 are 
given by 


*o = 

Fi = 
F 2 = 
F 3 = 


t [g-ltikx+ky+kz) 

-1- e’ 4 “ 

{kx — kx 

-k z ) 


+e ! 4 “( k * 

~\~ky — k z 

+ 

^ (~k x 

-k y +k z 

>]. 

4Aso sin 

(t) 

cos 

<*) 

— COS 

,ak z 

4Aso sin 

(t) 

cos 

{^-) 

— cos 

uk x 

4Aso sin 


cos 

dk x \ 

— cos 

■ Ciky 


)J. 

) 1 -( 121 ) 


This Hamiltonian exhibits 3D bulk Dirac points at three 
inequivalent X points X r = 2nf/a where r = a :,y,z. 
Each Dirac point at X r is protected by the 4-fold screw 
rotation about r axis. 

To understand the role of the screw rotation, let us de¬ 
scribe the symmetry of the system. The space group of 
the diamond lattice is FdAm{p 7 h ), which contains the 24 
symmorphic elements of tetrahedral point group 43 m(Td) 
and 24 non-symmorphic elements. The non-symmorphic 
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elements are obtained by compounding each symmor- 
phic symmetry operation of Td with a translation along 
td. = (f, f, f), which takes one sublattice site to another 
inequivalent sublattice site. On the other hand, the sym- 
morphic symmetry operation connects sites belonging to 
the same sublattice. 

At a generic point A = (0,0, k z ) on the k z axis, the 
system has C^ v symmetry that is composed of 5 different 
symmetry classes with elements {-E7|0}, {C\ |0}, 2 {C4\td}, 
2{iC\\td}, 2{iC' 2 \ti}, respectively. It is worth to note that 
some symmetry elements contain a partial lattice transla¬ 
tion td = (f, f, f) which is the characteristic property of 
the non-symmorphic nature of the diamond lattice space 
group. This contrasts with the case of conventional sym- 
morphic cubic lattices such as simple cubic (sc), face cen¬ 
tered cubic (fee), body centered cubic (bcc) lattices where 
the system on the k z axis has the ordinary Cd v group 
containing only symmorphic point group operations^. 

By considering the symmetry of the Hamiltonian, one 
can easily find the matrix representation of the screw ro¬ 
tation Cip = {Ci z \td} = exp (-ik-td)T x exp(ij<r z ) where 
exp(— ik ■ td) represents the translation td of the Bloch 
state with the momentum k, t x indicates the sublattice 
change due to partial translation, and exp(*^cr z ) repre¬ 
sents the 4-fold rotation. Then it is straightforward to 
confirm that 

C A ,iH(k x , ky, k z )C^j = H(-k y , k x , k z ). (122) 

Since 2,3 = ImF 0 = 0 on the k z axis, the Hamilto¬ 
nian on the k z axis becomes 

H(k z ) = 4£cos (123) 

from which we obtain one eigenvalue E+(k) = 

4tcos(|afcz) with the corresponding eigenvectors 


1 

f 1 ) 

0 

, x 1 

(o\ 

1 

71 

1 

■ lM = 7i 

0 


W 


W 


and the other eigenvalue E_(k) = —4tcos(jafc z ) with 
the eigenvectors 



( 1 ^ 


( 0 \ 


0 

-1 

• '^> = 71 

1 

0 


V 0 ) 


V 1 / 


Let us note that IV'+i)) ^+ 2 ), IV'-i); ^- 2 ) are also the 
eigenstates of 64,1 with the eigenvalues exp(—+ ij), 
exp(—- if), exp(-*^ + if), exp (—if- - if), 
respectively. Since the system is f periodic along the 
k z direction, a single Dirac point protected by C^i 
can be created at the Brillouin zone boundary. It is 
straightforward to see that these four eigenvalues have 


the same form as Eq. hence satisfy the condition 

for the band crossing at the Brillouin zone boundary with 
k z = ±27r/a. Also the band structure of the system along 
the k z axis is consistent with Fig. [ 8 ] 

Before we close this subsection, let us briefly perform 
a group theoretical analysis at the X point. The lit- 

_ x 

tie co-group G at the X point with the momentum 
k = (0,0, f) is generated by three symmetry operators, 

i.e., the four-fold screw rotation C^i about the z axis, 
the inversion P, the two-fold rotation about the x axis 
C 2 a) 27 ’ 43 - One interesting property of the diamond lat¬ 
tice is that the inversion symmetry P also accompanies 
a partial translation td- Thus, it is more suitable to use 
the notation P = {P\td} to indicate the partial trans¬ 
lation associated with the inversion. Here we describe 
two interesting physical consequences resulting from the 
non-symmorphic nature of P. 

Firstly, the partial translation involved in P does not 
affect the^commutation relation between P and a screw 
rotation Cpi,q on the rotation axis. This fact can be 
easily understood by considering the coordinate tranfor- 
mation under the combination of Cn j(? = {C/v|T g } and 
P = {P|£p}. Assuming the z axis is the screw rotation 
axis, PCN,q transforms the coordinate ( x,y,z ) to 

( % ^q,x T tp^x-. y T Z “b tp, z ), 

(124) 

where x' and y' are rotated coordinates satisfying x' + 
iy' = {x + iy) exp(*27r/A). Thus we obtain 

PC N ,q = {PC N I -Tq + tp}. (125) 

On the other hand, CN,qP transforms (x,y,z) to 

( — X ' + T q,x + t'p x - —y' + T q,y + tpyl ~ Z + T q,z + tp,z)i 

(126) 

where t' Px + it' Py = (tp tX + itp tV ) exp(*27r/A r ). Thus we 
obtain 

C N ,qP = {PC N \Tq + t'p}. (127) 

Now let us consider a Bloch state \k z ) on the rotation 
axis with the momentum k = (0, 0, k z ). We find 

PC N}q \k z ) = exp [ik z (-T q , z + t PtZ )} PC N \k z ), 
C N , q P\k z ) = exp [ ik z (T q , z + t P>z )] PC N \k z ), (128) 

which shows that the partial translation tp tZ associated 
with the inversion only provides an overall phase factor, 
and does not affect the commutation relation whereas r qtZ 
does. Therefore our theory can also be applied to systems 
with the inversion P accompanying a partial translation, 
as long as the Dirac point is located on the rotation axis. 

However, when the Dirac point is located away from 
the rotation axis, the translation t P can cause nontrivial 
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physical consequence as well. For instance, in the k z = 0 
plane, perpendicular to the rotation axis, a Bloch state 
\k x ,k y ) satisfies, 

PCN,q\k xl ky } 

= exp [ - ik' x (Tq iX - t P>x ) - ik' y (T q ,y - t P ^)\PC N \k x , k y ), 

(129) 

and 

CN,qP\kxi hy) 

= exp \ik' x (T qjX +f Px ) + ik' y (T qt y +t' Py )\PC N \k x ,k y ). 

(130) 

where k' x + ik' y = ( k x +ik y ) exp(?'27r/./V) and we have used 
the relation k' x t' P x + k' y t' Py = k x tp^ x + k y tp tV . The point 
is that since (t PtX ,t Pj y ) ^ {t'p x ^'p y ) due to the rota¬ 
tion, the partial translation tp associated with the inver¬ 
sion can also modify the commutation relation between 
the inversion and the rotation symmetries. Because of 
this, even a symmorphic rotation symmetry, which is not 
accompanied by a translation, can create a Dirac point 
away from the rotation axis when it is combined with the 
non-symmorphic inversion symmetry. 

For illustration, let us consider the commutation re¬ 
lation between C 22 and P = {P\tp} where tp = td = 
(|, |, j). The Bloch state \k x ,k y ) satisfies, 

PC 2 Z I k X : ky ) 

= exp [ik x (-tp <x ) + iky(-tp )y )\PC 2z \k x ,k y ), (131) 
and 

G-2zP\k x ,ky) 

= exp [ik x (t PiX ) + ik y (tp ty )]PC 2 z\k x ,ky). (132) 

Hence at the high symmetry momentum such as k\ = 
(qf, 0,0) or k 2 = (0,^,0) where exp (ik x t P}X + 

ikytp, y ) = i, we obtain {P,C 2z } = 0. This anti¬ 
commutation relation can create a stable Dirac point at 
k\ and k 2 , which is again confirmed by iF-theory analysis 
in the Appendix IB 31 Therefore, although the symmor¬ 
phic C 2 symmetry cannot support a stable Dirac point on 
the rotation axis, it can create a stable Dirac point at the 
zone boundary in the plane perpendicular to the rotation 
axis, when the C 2 is combined with a non-symmorphic 
inversion symmetry P. 

In fact, as noted before, there are three symmetry gen¬ 
erators, P, 64 , 1 , C 2x , at the X point with k = (0, 0, 
in the diamond lattice. According to the previous discus¬ 
sion, the Dirac point at X can be protected not only by 
P and 64,1 satisfying {P, 64 , 1 } = 0 but also by P and 
C 2x , which are also anti-commuting {P, C 2x j = 0. Be¬ 
cause of the high crystalline symmetry, the Dirac point 
in the diamond lattice is protected by multiple pairs of 
symmetry operators^!. 


VI. DISCUSSION 

To sum up, we have studied the topological charge of 
3D Dirac semimetals protected by the time-reversal, the 
inversion, and the rotation symmetries. Consideration of 
topological charges naturally leads to two different classes 
of Dirac semimetals, which is consistent with the previ¬ 
ous observation based on the symmetry constrained min¬ 
imal Hamiltonian analysis^. Class I Dirac semimetals 
are protected by an ordinary symmorphic rotation sym¬ 
metry which commutes with the inversion. Since each 
eigenstate carries a quantized rotation eigenvalue on the 
rotation axis, Dirac points should form a pair having the 
opposite topological charges when the system is periodic 
along the rotation axis. On the other hand, class II Dirac 
semimetals are associated with non-symmorphic screw 
rotation symmetries. The eigenvalue of a screw rotation 
is not quantized on the rotation axis due to the phase fac¬ 
tor induced by a partial lattice translation, which enables 
to create a single isolated Dirac point at the Brillouin 
zone boundary. 

The nonzero topological charge of a Dirac point not 
only guarantees the stability of the gap-closing point, 
but can trigger new types of Lifshitz transitions. For in¬ 
stance, when two Dirac points merge at the same momen¬ 
tum, the topological charge of the merging point is given 
by the summation of their topological charges. Since the 
energy dispersion around the gap-closing point strongly 
depends on its topological charge, such a merging tran¬ 
sition can generate intriguing nodal quasi-particles with 
novel physical properties'^. Moreover, the presence 
of a quantized topological charge can be a source of 
new topological responses. For instance, it is well-known 
that the nonzero monopole charge of Weyl points induces 
novel topological responses in Weyl semimetalsii. Recent 
theoretical studies of interesting topological responses in 
Dirac semimetal s 48 ' 49 may imply nontrivial role of topo¬ 
logical charges in these systems. 

Up to now, two materials (NaaBi and CdaAs 2 ) belong¬ 
ing to the class I are discovered and extensively studied 
whereas class II Dirac semimetals are not uncovered yet. 
Though there are some hypothetical candidate materi¬ 
als proposed theoreticall y 25 ' 26 , all of them are chemically 
unstable because the metallic ion in each candidate com¬ 
pound is required to have a lone-pair valence electron to 
locate the Fermi level at the Dirac points. In this re¬ 
spect, synthesizing class II Dirac semimetals is a challeng¬ 
ing problem in material science which should be properly 
addressed in near future. 

We believe that class II Dirac semimetals are as impor¬ 
tant as class I Dirac semimetals in the following sense. In 
the case of class I Dirac semimetals, the Dirac points are 
created by a band inversion, hence the intrinsic proper¬ 
ties of 3D Dirac particles can be observed only within the 
narrow energy scale associated with the band inversion^. 
Because of this, if the competing energy scales, such as 
the Fermi energy due to doped carriers or the disorder- 
induced broadening, are comparable to the band inver- 
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sion energy, the intrinsic properties of 3D Dirac particles 
can be easily washed out. However, in the case of class II 
Dirac semimetals, the energy scale of the Dirac dispersion 
is simply given by the bandwidth of the system (roughly 
in the order of a few eV), which obviously provides a 
better playground to study the intrinsic properties of 3D 
Dirac particles. 

Secondly, we would like to draw attention to class II 
Dirac semimetals as potential novel topological states. 
The fact that the Fu-Kane-Mele model is a canonical 
model to construct a 3D Z 2 topological insulator implies 
the intrinsic topological nature of the associated Dirac 
semimetal state. In particular, in the present paper, we 
have demonstrated that the presence of a single Dirac 
point on the rotation axis is unnatural in consideration 
of the Nielsen-Ninomiya theorem, and, in fact, the pro¬ 
jective nature of the screw rotation symmetry plays an 
essential role to circumvent the doubling of Dirac points. 
Although the discrete nature of the rotation symme¬ 
try should be distinct from the continuous U(l) symme¬ 
try associated with the original Nielsen-Ninomiya theo¬ 
rem, the mechanism leading to circumventing the fermion 
number doubling shares the common origin, i.e., assign¬ 
ing a non-quantized quantum number to fermions. To re¬ 
veal the topological properties of class II Dirac semimet¬ 
als would definitely be an exciting theoretical problem 
which we leave for future studies. 

Finally, we would like to point out that there are a class 
of Dirac semimetals which are not completely treated in 
our classification scheme. Let us note that, in both class 
I and II Dirac semimetals considered in the present work, 
Dirac points are located on the rotation axis. However, 
rotation symmetries can also protect a Dirac point which 
is away from the rotation axis. For instance, we have 
shown in Sec. IV FI that the symmorphic C 2 rotation can 
protect a Dirac point which is not on the rotation axis, 
when it is combined with a non-symmorphic inversion 
symmetry. Moreover, the tight-binding model on a hep 
lattice considered in Sec. IV El indicates that glide mir¬ 
ror symmetries can also give rise to symmetry-protected 
Dirac points. To find a systematic way to classify these 
different types of Dirac semimetals would also be an im¬ 
portant problem for future research. 


ACKNOWLEDGMENTS 

BJY is supported from the Japan Society for the Pro¬ 
motion of Science (JSPS) through the ‘Funding Program 
for World-Leading Innovative R&D on Science and Tech¬ 
nology (FIRST Program), and Grant-in-Aids for Scien¬ 
tific Research (Kiban (S), No. 24224009) from the Min¬ 
istry of Education, Culture, Sports, Science and Technol¬ 
ogy (MEXT). AF is grateful for support by Grants-in-Aid 
from the Japan Society for Promotion of Science (Grant 
NO.15K05141) and by the RIKEN iTHES Project. 


Appendix A: Absence of a stable Dirac point in 
systems with time-reversal and inversion symmetries 

only 

Here we prove that stable Dirac points do not exist 
in systems having only the time-reversal (T) and inver¬ 
sion (P) symmetries. For this, we distinguish two cases: 
one is when the Dirac point is located at a generic mo¬ 
mentum point, and the other is when the Dirac point is 
located at a time-reversal invariant momentum (TRIM). 
This distinction is necessary because the symmetry as¬ 
sociated with the Dirac point differs depending on the 
position of the Dirac point in the momentum space. In 
each case, the stability of a Dirac point is determined by 
using Jf-theor y 50 i 51 . 

1. A Dirac point located at a generic momentum 

When the Dirac point locates at a generic momentum, 
the combination of T and P is the only symmetry satis¬ 
fied around the Dirac point. In general, a PT symmetric 
system satisfies 

(PT)H(k)(PT)~ 1 = H{k), (Al) 

where the anti-unitary PT symmetry satisfies (PT) 2 = 
— 1 , which is coming from 

P 2 = 1 , T 2 = -1, [P, T] = 0, (A2) 

in electronic systems. 

The stability of the Dirac point can be understood by 
using K theory approach. We consider a Dirac point lo¬ 
cating at a generic momentum k Q = (fc°,fc°,fc°). The 
effective Hamiltonian describing the low energy excita¬ 
tion around the Dirac point is given by 

H d = (k x - k°)'jx + ( k y - ky)"f y + (k z - k°) 7 * + my 0 , 

where 'jo,x,y,z are gamma matrices satisfying the an¬ 
ticommutation relations { 7 ^, 7 j} = 2 5ij, and m indi¬ 
cates a possible symmetry-preserving Dirac mass term. 
The presence (absence) of the symmetry-preserving Dirac 
mass term m indicates the instability (stability) of the 
Dirac point. From Eq. m, we obtain 

hx,PT\ = [7 v,PT] = [ 7z , PT] = [70 ,PT] = 0. (A3) 

To confirm the existence or absence of the Dirac mass 
term m, let us define a real Clifford algebra generated by 
PT and 70 , x ,y, z , which is given by 

Clo, 4 8 Cl 2 , 0 = {; lx,ly, 7 *> 7o} <8 {PT, JPT; },(A4) 

where C7 Pi9 indicates a real Clifford algebra with p + q 
generators {ei, e 2 ,..., e p ; e p+ i, e p+2 ,..., e p+q } satisfying 

{e i ,e j } = 0, i^j 

— 1 , 1 < * < p, 

+1, p+l<i<p + q. 
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The algebra in Eq. (I ATI) can be easily obtained by con- where the matrices 70 tX ,y,z satisfy { 74 , 7 , 7 } = 2 Sij. T and 
sidering the following relations: P symmetries require 


(0 [ 7 x,PT\ = [7 y,PT] = [7 z,PT] = [ 70 , PT] = 0, 

(ii) [ 7 a,, J] = [ 7 ^ , J] = [l z , J) = [ 70 , J] = 0, 

(m){7i,7 j} = 2 Sij, 

(iv) {PT, JPT} = 0, (PT) 2 = (JPT) 2 = - 1 , (A5) 

where the symbol J indicating the pure imaginary num¬ 
ber i is introduced to construct a real Clifford algebra. 

The existence or absence of the Dirac mass 77170 can be 
judged by considering the following extension problem: 

{;7x,7!/} ® {PT, JPT;} 

—>{; 7 x ,7 yilz}®{PT,JPT]}, (A 6 ) 


{lx , T} = { 7y , T} = { 7z , T} = [ 7o , T] = 0, 

{ 7 x; P} = { 7 ?y- P} = {lz,P} = [70, P] = 0 , (A 12 ) 

where 

P 2 = 1,T 2 = -l,[T,P] = 0. (A13) 

The Clifford algebra generated by jo,x, y , z> T, P, J is 
given by 

{T, JT, J 70 ; 7x, ly , 7* 1 Plxlylz }, (A14) 

The existence or absence of the Dirac mass 77170 can be 
judged by considering the following extension problem: 


i.e., 


C 7 0 ,2 ® < 722,0 ~t C/0,3 ® <7/2,0- (A 7 ) 


{ J7 1 ; 7x, ly , 7*. -P7x7y7z} 

-t {T, JT;7x,7j / ,7 z ,P7ai7y7 z } (A15) 


Namely, the topological classification of 7 Z determines 
the topological nature of the Dirac point. This is because 
the topologically trivial classification of y z implies the 
existence of another gamma matrix such as 70 which an¬ 
ticommutes with the three generators l x , y ,z (i.e., a mass 
term exists) whereas the topologically nontrivial classifi¬ 
cation of y z implies the absence of 70 , thus the topolog¬ 
ically nontrivial nature of the Dirac point. Generally, in 
the classification scheme with Clifford algebera, the exis¬ 
tence condition of a particular generator (Dirac mass 
term) is equivalent to the classification of another gen¬ 
erator of the same type as e* in the Clifford algebra in 
which ei is removed^ 

The extension problem 


<7/0,2 ® <7/2,0 ~ f <7/0,3 ® <7/2,0 

(AS) 

is equivalent to 


< 7 / 4,0 —)► <7/5,0 

(A9) 

due to the relation 


Clp, q G C/2,0 — <7/4+2, p- 

(A10) 

Since the classifying space for the extension 
< 7 / p +i ,4 is given by P p + 2-q, the classifying 

<7 /p , 9 — 7 
space for 


the extension C/ 4,0 —> < 7 / 5,0 is Re = Sp(n)/U(n) with 
a sufficiently large integer n. Since 7r 0 (P6) = 0, the 
space of possible representation for y z is singly connected. 
Namely, a Dirac mass term always exists, hence the Dirac 
point is unstable. 

2. A Dirac point locating at a TRIM 

On the other hand, when the Dirac point locates at 
a TRIM, both P and T are the symmetry of the Dirac 
point. To understand the stability of the Dirac point, we 
consider the following Dirac Hamiltonian 

Hd = k x i x + ky'Jy + fc z 7 z + 77770, (All) 


i.e., 

C/ 1,4 —^ < 7 / 2 , 4 , (A16) 

the corresponding classifying space is R p + 2 -q = R -1 — 
P7. Since 7To(P7 ) = 0 , a Dirac mass term always exists, 
hence the Dirac point is unstable. Therefore, indepen¬ 
dent of the location of the Dirac point in the momentum 
space, the system with only T and P symmetries cannot 
support a stable Dirac point. 

Appendix B: The stability of Dirac points in C 2 
symmetric systems 

In C 2 invariant systems, since all the symmetry opera¬ 
tors squared become ±1, the classification scheme based 
on Clifford algebras can be applied. 

1. When [C 2 , P] = 0 

In systems with the symmorphic C 2 rotation satisfying 
(C 2) 2 = — 1, the Dirac point can be located at a generic 
momentum fco = (0,0, fc°) on the rotation axis (z axis). 
We consider the following massive Dirac Hamiltonian, 

H(k) = k x ^ x + k y ly + (k z - k° z ) q z + 77170, (Bl) 

which satisfy the following relations, 

(PT)P(fe)(PT)- 1 = H{k), (B2) 

C 2 H(k x , ky.k^r 1 = H(-k x , -k y , k z ). (B3) 

From this, we obtain 

[7,, PT] = [ 7l „ PT] = [ 7z , PT] = [70, PT] = 0 , (B4) 

{7x, c 2 j = {ly, c 2 j = [ 7 z, <7 2 ] = [70, c 2 ] = 0 . (B5) 

Then we can determine the Clifford algebra generated 
by the gamma matrices in the Dirac Hamiltonian and 
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the relevant symmetry operators. The resulting Clifford 
algebra is 

C7 6 ,o 0 C7 0 ,i 

= {PT,JPT,J'y x ,Jj y ,Jj z ,Jjo;}0{;7x7yC 2 }. (B6) 
The relevant extension problem is 

Cl 4 ,o 0 Clo,i Cl B ,o 0 Clo,i, (B7) 

for which the classifying space is given by R e x R e . From 
its zeroth homotopy group, we find the topological charge 
as 7To(i?6 x P&) = 0. Therefore a Dirac point cannot carry 
a nontrivial topological charge, which is consistent with 
the absence of a topological invariant found before. (See 
Table I.) 

2. When {C 2 ,P} = 0 

Now we consider a two-fold screw rotation C 2 satisfy¬ 
ing (C 2) 2 = 1. Since the Dirac point locates at a TRIM, 
both the P and the T are the symmetry of the Dirac 
point. The commutation relations relevant to this prob¬ 
lem are as follows. 

{ 7 „ T} = { 7y , T} = { 7z , T} = [ 7o , T] = 0, 

{ 7 *, P} = { 7 y, P} = { 7 zi P} = [ 7 o ,P\ = 0, 

{7x, C 2 } = { 7 y> C 2 } = [7 zj C 2 ] = [ 70 , C 2 ] = 0, (B8) 

and 

[T,P] = 0, [C 2 ,T\ = 0, {C 2 ,P} = 0, (B9) 

in which 

P 2 = 1, T 2 = —1, (C 2) 2 = 1. (B10) 

The relevant Clifford algebra of gamma matrices and 
symmetry operators is given by 

Cl 3 ,5 = {T, JT , J7a-,7 x ,7y,7 z ,P7 x 7y7z,C 2 P7z}- 

(Bll) 

The existence condition of the Dirac mass term is deter¬ 
mined by the extension problem 

Cll,5 “A C^2,5; (B12) 

for which the classifying space is Rq. From tt 0 (R 6 ) = 0, 
we see that the C 2 symmetry cannot protect a Dirac point 
consistent with Table III. 

3. When {C 2 , P} = 0 

This is relevant to the case when the Dirac point is 
located at a TRIM in the plane perpendicular to the 


rotation axis. The commutation relations relevant to this 
problem are as follows. 

(7„ T} = {■ ly , T} = { 7z , T} = [70, T} = 0, 

{7x, P} = {7y: P} = {7z, P} = [70, P\ = 0, 

{7x, C 2 } = { 7 y , C 2 } = [ 7 z,C 2 ] = [ 7 o, C 2 ] = 0, (B13) 

and 

[T,P] = 0, [C 2 ,T] = 0, {C 2 ,P} = 0, (B14) 

in which 

p2 = _ 1; r 2 = _ 1; (C 2 )2 = - 1 . (B15) 

Then the relevant Clifford algebra of gamma matrices 
and symmetry operators is given by 

CI 3 ,5 = {T, JT , J70iP7x7ylz\7x,7y,lziC 2 P7 z }. 

(B16) 

The existence condition of the Dirac mass term is deter¬ 
mined by the extension problem 

Cl 2 ,4 i d 3l 4i (B17) 

for which the classifying space is Rq. From 7To(i?o) = 2, 
we see that the symmorphic C 2 symmetry can pro¬ 
tect a Dirac point when it is combined with the non- 
symmorphic inversion symmetry P. Here the location 
of the Dirac point is not on the rotation axis but at a 
TRIM on the plane perpendicular to the rotation axis 
because the anti-commutation relation {C 2 ,P} = 0 can 
be satisfied only away from the rotation axis. 

Appendix C: The stability of 2D Dirac points in 
systems with two-fold screw rotations 

Recently, Young and Kane proposed a theory^ about 
2D Dirac points located at a TRIM on the Brilluoin zone 
boundary. One interesting finding in their work is that 
a two-folcl screw rotation can protect a Dirac point at a 
TRIM on the rotation axis, which is forbidden in the case 
of 3D systems. To confirm their claim, let us check the 
stability of the Dirac point by using K theory. 

Let us consider a Dirac Hamiltonian at the zone bound¬ 
ary. 

H(k) = k x 7 x + k v 7 v + m 7 0 , (Cl) 

where 7o,x,y are mutually anti-commuting. Under P, T, 
C 2y symmetry satisfying 

[P,T\ = [C 2y ,T] = 0, {C 2yi P} = 0, 

P 2 = 1, T 2 = -1, Cly = 1, (C2) 

the gamma matrices satisfy 

{7 x,T} = { 7 y, T} = [ 70 , T) = 0, 

{7*, P} = {ly, P} = [7o, P] = 0, 

{lx , C 2 y} = [7y, C 2 y] = [ 7 0, C 2 y\ = 0, (C3) 
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The relevant Clifford algebra of gamma matrices and 
symmetry operators is given by 

Cl4,2 ® Cl\ o — {T, TJ,J yo> C2y r )x i Tx 5 Ty } ® { C^fx^fy > } 

(C4) 

The existence of the Dirac mass term myo is determined 
by the extension problem 

C? 2,2 ® Clifi —t Cls t 2 ® Clifi. (C5) 


Since the extra generator P^jxfy commutes with all the 
other generator, and satisfies {P^xly) 2 = — 1 , the above 
extension problem is rearranged in the following way, 

Ch -t Cl 6 , (C 6 ) 

for which the classifying space is C 4 = Cq. Namely, the 
extra generator P^xAy deforms the original real Clifford 
algebra extension problem to a complex Clifford algebra 
extension problem^. From 7 To(Co) = Z, we see that C 2 
can protect a 2D Dirac point with topological charge Z 
on the rotation axis. 
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